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Ch ■ Abstract 

Y \ We study via mean- field theory the pairing between Fermions of different masses, especially at 

Qh' the unitary limit. At equal populations, the thermodynamic properties are identical with the equal 

mass case provided an appropriate rescaling is made. At unequal populations, for sufficiently light 
majority species, the system does not phase separate. For sufficiently heavy majority species, the 
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^^ ' phase separated normal phase have a density larger than that of the superfluid. For atoms in 

O ' harmonic traps, the density profiles for unequal mass Fermions can be drastically different from 
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their equal-mass counterparts. 
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Cooper pairing between Fermions due to the presence of an attractive interaction is the 
essential physics for most, if not all, superconductors. In recent years, Feshbach resonances 
[l| provide us with a wonderful opportunity to further study strongly interacting Fermi 
systems. For wide resonances, the Feshbach resonance between the atoms can be viewed as 
providing an attractive pairing interaction, with the pairing strength increases from positive 
detuning towards negative detuning. For the case of equal populations of two Fermionic 
species and an s-wave interaction between them, it was predicted long time ago j^ that 
the system can smoothly crossover from a BCS superfiuid to a Bose Einstein condensate of 
tightly bound pairs. Recent experiments have provided a confirmation of this picture, and 
many thermodynamical and dynamical properties have been measured. |3| There is also an 
explosion of theoretical activities studying this system, j^ 

In the case where the populations of the Fermions are unequal, the smooth crossover 
discussed above is known to be destroyed J5| . Recent experiments p have shown that phase 
separation occurs near resonance, so that the superfiuid appears near the center of the trap 
and a normal fluid of unequal population density appears outside the superfiuid core. [3] 

In the literature cited above, the Fermi atoms involved are identical except their internal 
electronic states. Hence they have identical masses. However, Feshbach resonances are not 
confined to the same atom species ^Sj. This thus raises the possibility of studying Fermion 
pairing between two atomic species of very different masses, with either equal or unequal 
populations. This problem can also have important implications in other areas of physics 
such as quark matter [2| and superconductivity in general. 

In this paper, we thus study two Fermionic species of unequal masses with an interspecies 



short-ranged s-wave attractive interaction. First we consider t 



le equal population case. We 







show that, at zero temperature and within mean-field theory 2|, the thermodynamic prop- 
erties of this system have a one-to-one correspondence with those of the equal mass case, 
provided we rescale the energies appropriately. Next we consider the case of unequal pop- 
ulation, but confining ourselves to the unitary regime where the scattering length diverges. 
Defining 7 as the ratio of the majority atom mass to the minority mass, we show that the 
system phase separates for 7 > 7^;/ ^ 0.13 but remains homogeneous for 7 < '-y^. For suffi- 
ciently light majority particles, it is energetically costly to phase separate out the majority 
particles since a high Fermi energy of the majority species will result. For 7 > 7//, we show 
that the system phase separates into a superfiuid with equal population density and a nor- 



inal fluid with unequal population density. However, there is dramatic difference between 
the cases where 7 is larger or less than a value 7* ~ 3.9. For 7h < 7 < 7*, the density of the 
normal fluid is less than the superfluid, as one may expect due to the attractive interaction. 
However, for 7 > 7*, the normal fluid is actually denser than the superfluid. We explain 
this by the larger density of states for the heavy majority component. We illustrate this 
difference by calculating the density proflle of an unequal population gas in an harmonic 
trap. 

We shall then consider at zero temperature two Fermionic species with a short-ranged 
interaction characterized by the s-wave scattering length a. 1/a increases from —00 to +cx3 
with the strength of the pairing interaction. We shall call the two species up and down, 
with masses 771^^, a =] or [. The effective mean-fleld Hamiltonian is -ffm/ ~ /i-f^i — /i^A^j, = 
^A:,.e.(A:)4/,>-E£[A*c_,^a^M+c-c]-[2^-Eg^]|Ap. Heree.(fc) = h^ky2m^-^^^ 
are the quasi-particle excitation energies for normal fermions relative to their respective 
chemical potentials /lo-, Cj:J^ are the annihilation operators for the Fermion a at momentum 
k, c.c. denotes the complex conjugate, rur is the reduced mass (m~^ = m^^ + mj^). A, the 
off-diagonal pairing potential, is determined by minimization of the expectation value of this 
Hamiltonian {i.e. the free energy [lO|) with respect to A. This gives, for equal populations, 
(see, e.g |ll|) 
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E„{k), the quasiparticle energies at wavevector k, are given by 



(1) 



EM ^ m^s^A^ + j(m±LA^^ + A. 
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where and — t— i- Deflning /i = (yU| + /i|)/2 and h = (/i| — yU|)/2 and noticing that 
^„{k) + i-M = 1^ - 2/x, one obtains 



E^{k) + E^{k) = 2^ 
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-/i)2 + |A|2 = 2E(A;) 
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The number equation is 
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where v{k), the usual coherence factor, is given by JjJ v'^{k) = J.]^ = p;\e ^^"^ hence 
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v\k) = -\l-^£l^LZ^\ (5) 

The equations (Hj) and (JD) determine A and /i for a given density. Using eqs Q and 
(jSj), we see that these equations are identical to those [2] of equal masses m^ = mi = m 
provided we substitute m -^ 2mr. Thus their solutions are identical after this replacement. 
(The solutions are thus independent of h so long as all quasiparticle energies Erj{k) > 0: 
This condition also guarantees that n^ = n^, see below). Therefore their solutions can be 
expressed as 



where t 

masses 



le dimensionless functions Ci,2 are identical to the corresponding ones for equal 
^. Here epr = kp/2mr and kp = (STr^n)^/^, where n is the total density. We have 
introduced a factor of 1/2 in eqs © and ((7j) since eprf^ = k\/{2m) if mg- = m. It follows 
immediately that other thermodynamical properties can be obtained similarly. For example, 
the energy density Es/V can be written as Es/V = n^Cz ( -^-^ ) where Cs is again the same 
function as the equal mass case. 

In particular, at unitarity (a = oo), we have yU = ^Ci(O) where Ci(0) ~ 0-59 in the 
present theory (and thus /i > 0). The magnitude of the order parameter A = ^C2(0) with 
(^2(0) ~ 0.68 and hence A/yU ~ 1.16. Using the fact that the only length scale at unitarity 
is kp^, (3(0) = |Ci(0) and thus Eg/V = |n^Ci(0)- At chemical potential /x, from eq ©, 
kp = [4mr/i/^i(0)] , hence the density is given by n = -^ [4mr/i/Ci(0)] . Since the free 
energy density is fls/V = Es/V — fin = —^rifi, therefore 

^ ^ 2_ (4rrvfV^ 

V 157r2 ^3/2 ^Q^ • i»J 

Mean field theory is known to produce a larger Ci(0) than Monte-Carlo and experiments 
(^ 0.44) for the equal mass case ^\^. It is however of interest to see to what extent Ci etc 
will depend on 7 in the unequal mass case in more accurate calculations. 



Next we consider unequal populations but confine ourselves to unitarity (1/a = 0). We 
shall call t (i) the majority (minority) species. The self-consistent order parameter equation 
reads 
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(9) 



E^{k) + Ei{k) h^k^ 

Due to the appearance of the Fermi functions / with arguments Erj{k), the solution to this 
equation does not depend on the reduced mass alone, ( nt > n^ requires that some E^{k)^s 
are negative, since n^ — rii = y^kf{E^), see, e.g. p, |ll|). By rescaling all energies to 
/i, it is easy to see that eq Q gives A//i as a function of /i//i for each 7 = m^/m>^. The 
non-trivial (A 7^ 0) solutions to eq (jH)) consists of two branches (see inset of FigQ]). One 
branch has A//i ^1.16 independent of /i//i and 7. This branch corresponds to E(j{k) > 
for all k and a, i.e., no quasiparticles (and hence equal populations). This solution has 
already been described in the last paragraph above. Another branch, corresponding to 
^d = n^ — ^i y^ 0, gives a /i//i and 7 dependent A//x. For 7 < 0.13, A//i decreases with 
increasing h/fj.. We have verified that this branch is a free energy [lO| minimum, and thus 
represents a stable homogeneous superfluid. (We have also verified that this state has only 



one Fermi surface, and hence it does not correspond to a breached-pair state 
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However, 



for 7 > 0.13, A/yU increases with increasing h/fi (and correspondingly, for some range o 
A//i is multi- valued). This branch is a generalization of the solution obtained by Sarma 12 1 
in the BCS limit for 7 = 1. We found that this solution corresponds to a free energy 10 1 
relative maximum and is thus unstable. 

For a homogeneous system with 7 > 0.13, the likely scenario is that the system phase 
separates, as |6^ in the case of 7 = 1. (We ignore the more exotic possibilities jlJ] here). To 
characterize this phase, we seek the point where the free energy density of the completely 
paired superfluid state is equal to that of the normal fluid, i.e., ^^(/i, h)/V = fi7v(yU, h)/V . 
^s/'^ is independent of h and was already given before. Within our mean field theory, 
VLn /V is just the free energy density of an non-interacting Fermi gas, thus VLn /V = —j^ 
[(2mt)^/^(yu + hf/^ + (2m|)3/2(/i - hf/^H] where 7^ = 1 if /x - /i > and = otherwise. 
Equating the free energies, we obtain an equation for h/fi. Denoting this value as {h/fj,)c, 



we get (using /i > at resonance) 
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,/^/c/ V ^ / V v^. 

For h/fi < {h/fi)c, we have a paired superfluid, whereas for h/fi > {h/fi)c, we have a 
normal fluid. At h/fi = {h/fi)c, the system phase separates into a superfluid region and a 
normal region. The dependence of {h/fi)c on 7 is given in the main part of Fig^ {h/^)c > 1 
for small 7, decreases with increasing 7, equals unity at 7 = 7p = Cr^(O) — 1 ~ 0.69, changes 
sign at 7 = 7,, ~ 3.9, and is negative for 7 > 7,,. (One should not be alarmed by this sign 
change: the concentration of up particles in the normal state is larger than that of the down 
particles if m-f{fi + h) > mj(/i — h), i.e., - > j^-) This behavior is the result of unequal 
masses of the two species. This is particularly clear if we examine the behavior of {h/fi)c for 
large 7. In this limit, the solution to eq (fTUj) is given asymptotically by {h/ ix)^ = — l + ci/7^/^ 
where Ci = 2 [(Ci(0))^^^^ — l] • The dependence on 7 is due to the different density of 
states between the two species. With Ci(0) ~ 0.59, we find Ci ~ 2.16. This asymptotic is 
plotted also in Fig ^ which gives a reasonable description to the actual solution. 

The ratio of the minority to majority density is {njn^)^ = 

H [(1 - (/i//i)c)/(7(l + {h/ii)c))f^'^. {nd/n)N can be obtained via (1 - {nJni^)N)/{l + 
{njn^)]^). This ratio as a function of 7 is plotted in Fig |21 {nd/n)iy = 1 for 7 < ■jp. 
There the only stable normal phase is the completely polarized phase. Since n^ = n-^ in 
the superfluid phase, {nd/n)^ is also the maximum allowed {nd/n) for the phase separated 
system. For < nd/n < {nd/n)^, the system phase separates into a superfluid and a normal 
state, whereas for {nd/n) > {nd/n)^ the system is in the normal state. These two different 
regions are labelled as PS and N in Fig|21 The dependence of {nd/n)N on 7 is a result of 
competition between the 7 dependence of {h/fi)c and that of the density of states. For large 
7, using the asymptotic form for {h/jj?)c discussed before, we get {njni)^ ~ (2/ci)^/^/7^/^, 
and thus {nd/n)^ ^ 1 — 2^/^/(c;^ 7'^/^). This asymptotic dependence is also plotted in Fig 
El The dependence of {ni/n'^)]\f on 7 is a combined result of (1 + {h/ fi)c) oc 1/7^/^ and the 
ratio of density of states of the down to up particles (oc 7~^). This competition also results 
in the non-monotonic behavior of {nd/n)^ on 7, with a minimum occuring at 7 slightly 
larger than 7^,. 

The density of the normal fluid is given by n^- = 

g^ [(2m|(/i + /i))^/^ + (2m|(/i — /i))'^/^7Y] and thus its ratio to that of the superfluid 



state is 
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The dependence of n^/ns as a function of 7 is as shown in Fig El n^/ns increases with 
increasing 7. For 7 < 7* ~ 3.9, UN/ns < 1, but for 7 > 7*, un/us > 1. At 7* where 
{h/ ji)c = 0, n^ /ns is exactly unity as can be verified easily by comparing eqs (fTUI) and (fTT|). 
For large 7, we have n^/ns ~ Ci (0)[ci 7'^''^ + 2'^/^]/2^/^, hence n^/ns increases with 7 
roughly as 7^/^. This dependence is again a combination of 7 dependence of {h/ix)c and the 
density of states factors. The asymptotic behavior is also plotted in Fig El 

The ratio nj^/ns and {na/n)^ above are directly reflected in the density profile of a 
trapped gas. Let us consider for simplicity an isotropic harmonic oscillator trap, with trap 
potential energy for the two species ^a^r"^. The local chemical potentials are given by 
/io-(r) = /i° — |ao-r^, and hence fi{r) = /i° — ^{a^ + ajr^ and h{r) = h° — ^{a-^ — ai)r'^. 
Generally h{r)/^{r) increases from the center towards outside. Therefore a superfluid core 
exists for r < Tc where r^ is defined by h{rc)/fi{rc) = {h/fi)c- In this region, the local density 
n{r) is given by (4:mrfi{r)Y''^ / (^rr'^Ci (0)) with n'^{r) = ni{r) = n{r)/2. For r > Tc we 
have a normal fluid. There na{r) = [2mo-(yu(r) ± /i(r))]'^/^/(67r^) provided fi{r) ± h{r) > 
and vanishes otherwise. As an example, we show in Fig JH the density profile for the case 
«! = «! = a, 7 = 6.7 (~ m(^°K)/?7i(^Li)). This figure clearly shows that the density 
at r = r+ is larger than that of r~ (with ratio ~ 1.2, c.f. Fig jSl also {ni/ni^)\j.+ ~ 0.12 
corresponding to {nd/n)]^ ~ 0.78, c.f. Fig El). For comparison, the density profile for the 
case of equal mass (7 = 1) is shown in the inset of Fig 01 which shows a density drop from 
the superfluid core to the normal region. We also note that in our present case, the minority 
particles have a larger cloud radius than the majority. This is because of the lighter mass 
of the minority species, so that in fact yU° > /i° and /i° < 0. 

In conclusion, we have studied the superfluid pairing between unequal mass Fermions, in 
particular at resonance. We show that, for equal populations, the thermodynamic properties 
of the system are the same as the equal mass case except for a simple rescaling, whereas for 
unequal populations, they depend crucially on the mass ratio of the two species. There is an 
abundance of stable or long-lived isotopes with very different masses even within the alkalis, 
ranging from ^H to ^^^Cs, and thus ample opportunities to test the present predications. For 
example, at resonance, a ^H(majority)- ^°K(niinority) mixture would not phase separate, 



whereas a ^'^K(niajority)- ^H (minority) mixture would. For ^Li(majority)-'^°K(minority) 
mixture, the system would phase separate with a density drop from the superfluid core to 
a (completely polarized) normal region, (in contrast to equal mass case, where the normal 
region is only partially polarized), whereas '^°K(majority)-^Li (minority) mixture would have 
a density jump. 

Very recently, another preprint [ij] appears which also deals with unequal mass Fermions. 
This research was supported by the NSC of Taiwan under grant numbers NSC94-2112-M- 
194-008 (STW), NSC94-2112-M-194-001 (CHP), and NSC94-2112-M-001-002 (SKY). 
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FIG. 1: (color online) (h/fi)^ as a function of 7 for 7 > 7/^. The dashed line is the asymptotic 
formula in text. Inset: A/fj, as a function of h/fi for 7 = 0.14 and 7 = 0.12. 
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FIG. 2: (color online) The bulk phase diagram at resonance. For rid/n = 0, the system is a super- 
fluid, whereas it is normal for Ud/n = 1. For < n^/n < 1, the phases are: H - the homogeneous 
phase, PS - phase separated state, N - normal phase. The dashed line is the asymptotic formula 
discussed in text. 
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FIG. 3: (color online) The ratio of n^ to ng at phase equilibrium for the phase separated region. 
The dashed line is the asymptotic formula discussed in text. 
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FIG. 4: (color online) Density profile at resonance for iV = 4.1 x 10^, N^ = 6.6 x 10'', N^/N = 0.16. 
«! = a| = Q, 7 = 6.7. Full line: total density, dots and dashed: n^^n) in the normal fluid region. 
Here r is normalized to / = l/[a{2mr)Y''^ ^"^^ ^^^ density n to /^^. Inset: density profile for the 
same N and N^ for 7 = 1. n\^[r) for the normal phase is very small and thus only barely visible. 
The density profiles with the same N^/N but A'^ncw = P]^ can be obtained by the scaling relation 
nnew(/3'/'r)=/?i/2^(r). 
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